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Abstract 

Aims. The Crab nebula displayed a large 7-ray flare on September 18, 2010. To more closely understand the origin of this 
phenomenon, we analyze the INTEGRAL (20-500 keV) and FERMI (0.1-300 GeV) data collected almost simultaneously 
during the flare. 

Methods. We divide the available data into three different sets, corresponding to the pre-flare period, the flare, and 
the subsequent quiescence. For each period, we perform timing and spectral analyses to differentiate between the 
contributions of the pulsar and from the surrounding nebula to the 7-ray luminosity. 

Results. No significant variations in the pulse profile and spectral characteristics are detected in the hard X-ray domain. 
In contrast, we identify three separate enhancements in the 7-ray flux lasting for about 12 hours and separated by an 
interval of about two days from each other. The spectral analysis shows that the flux enhancement, confined below 
~1 GeV, can be modelled by a power-law with a high energy exponential cut-off, where either the cut-off energy or the 
model normalization increased by a factor of ~ 5 relative to the pre-fiare emission. We also confirm that the 7-ray fiare 
is not pulsed. 

Conclusions. The timing and spectral analysis indicate that the 7-ray fiare is due to synchrotron emission from a very 
compact Pevatron located in the region of interaction between the pulsar wind and the surrounding nebula. These are 
the highest electron energies ever measured in a cosmic accelerator. The spectral properties of the flare are interpreted 
in the framework of a relativistically moving emitter and/or a harder emitting electron population. 

Key words. Acceleration of particles - Astroparticle physics - Magnetic fields - (Stars;) pulsars: individual: Crab - 
Gamma rays: stars - X-rays: stars 



1. Introduction 

With an integrated luminosity of about 5 x 10'^* ergs^^ 



and 



a distance of ~2 kpc, the Crab supernova remnant is very 
bright from the radio domain to TeV energies (see e.g., 
for a review). 



Hester 2008 



It is powered by a pulsar spin- 
ning on its axis in about 33 ms that injects energetic elec- 
trons into the surrounding nebula. Nearly all the nebular 
emission up to 0.4 GeV is believed to be produced by syn- 
chrotron cooling of these electrons in an average magnetic 
field of ^ 300 /iG. At higher energies, inverse Compton (IC) 
cooling dominates. 

The integrated high-energy flux of the nebula and the 
pulsar has been remarkably stable over the past few decades 
and the object is indeed used as a c alibration source in sev- 
eral experiments (but see j Wilson-Hodge et al. 2010 for the 
first report of a secular X-ray trend). On 2010 September 22 



the AGILE collaboration (Tavani et al. 112010 1 reported the 



first 7-ray flare from a source positionally consistent with 
the Crab, during which the flux above 100 MeV was nearly 
double its normal value. The 7— ray flare from the direc- 
tion of the Crab was confirmed by the Fermi collaboration 
( Buehler et aL][20Tol ). 

During a period partially covering the 7-ray flare, the 
Crab region was observed by the INTEGRAL satellite and 
the Swift/BAT telescope during its routine sky survey. No 



statistically signiflcant increase in the Crab flux was ob- 
served in the IBIS/ISGRI light curves between 20 and 400 
keV, as well as in the Swift/BAT 15 - 50keV flux at a 



level of 5% at the la confidence level ( Ferrigno et al.pOlO 
Markwardt et al.||2010[ ). 



Swift /XRT observed the Crab region for Iks on 2010 
September 22 at 16:42 UT and did not reveal any signifi- 
cant _va£iatiOT_intl^^ spectrum, and pulse pro- 
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file (Evangelista et al. 2010). The ultraviolet and soft X- 
ray images excluded the presence of any bright unknown 
field object th at could have contributed to the 7-ray flux 
(Heinke 2010). The near-i nfrared Crab flux in the J and H 
bands was also constant ( jKanbach et al.pOlO l. Dedicated 
pointing performed by RXTE did not show any signiflcant 
chang e in the overall spectral p roperties in the 3-20 keV 
band ( Shaposhnikov et ar]|2010 ). 

A 5ks TOO Chandra observation was performed on 
2010 September 28 to monitor the morphology of the in- 
ner nebula. There are no particular variations with respect 
to the over 35 previous observations, with the possible ex- 
ception of an anomalous extension of a bright knot closer 
(down to 3") to the pulsa.r, als o noticed in one archival ob- 
servation ( Tennant et al.|[2010 ). A subsequent HST optical 
observation of the Crab conflrmed an increase in the emis- 
sion about 3 arcsec sout h-east of the pulsar with respect 
to archival observations ( Caraveo et al.pOlO ). However, it 
remains unclear whether this feature is related to the 7-ray 
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event. The Crab 7-ray flux returned to its usual level on 
2010 September 23, les s than a week after the onset of the 
flare ( |Hays et al.||2010[ ). 



ISGRI 20-40 keV 



2. Results 

2.1. INTEGRAL 

The INTEGRAL satellite is equipped with several high- 
energy instruments covering the energy range from a few 
keV to a few MeV. Her e, we concentrate on th e results from 
the spectrometer SPI (Ve drenne et al. 2003) and t he hard 
X-ray imag er IBIS/ISGRI ( [Ubertini et al.n2003l [Lebrun 
et al.|2003 1 . We divided the available observations into two 
periods: the pre-flare from 2010 September 12 10:32 UT to 
2010 September 18 05:56 UT (exposures: ISGRI 280 ks, SPI 
336 ks), and the flare from 2010 September 18 11:24 UT to 
2010 September 23 10:20 UT (exposures: ISGRI 89 ks, SPI 
109 ks). The gap between the two data sets is due to the 
uncertainty in the exact time of the flare onset. The data 
are analysed using the off-line science analysis so ftware ver- 
sion 9 provided by the ISDC (jCourvoisier et al.||2003). The 



data analysis is performed using the SPIROS package and 
the GEDSAT template background model, based on the count 
rate of the events saturating the Germanium detector. The 
spectra obtained for both data sets are equal within the 
errors. For the pre-flare set, a best-fit to the data is ob- 
tained using a bro ken power-law model with a fixed energy 
break at 100 keV ( [Jourdain fc Roques||2009[ ). The spectral 
indices are Fi = 2.142 ± 0.019 and "'2 = 2.194 ± 0.021 
(X^ = 42/47). For the data-set corresponding to the flare, 
a broken power-law model is also adequate, whose best-fit 
indices are Fi = 2.103 ± 0.038 and Fa = 2.149 ± 0.041 (x^ 
= 52/47). The total flux between 50 and 1000 keV is, re- 
spectively, 2.79tQ;°2 ^ srg cm~^ s^^ before the flare 
and 2.8llo'2e x 10~® erg cm~^ s~^ during the flare, which 
are equal within the uncertainties. 

While the absolute count rate in IBIS/ISGRI is subject 
to a few percent variation due to instrumental systematic 
effects, the background-subtracted pulse profile normalized 
to its average count rate can be considered a trustful esti- 
mator of the pulse shape. We can thus compare the relative 
contribution of the pulsar and nebular emissions before and 
during the 7-ray f lare. We exploited the metho d of [Segreto 
fc Ferrig no (2007) and used the ephemeris of Lyne et al. 
to extract background-subtracted pulse profiles, ac- 



(T993 



cumulated in 100 phase bins in the energy ranges 20-40 keV 
and 40-80 keV, and in 25 bins in the 80-150 and 150-500 keV 
energy ranges. 

The sum of the squared differences between the pulse 
profiles' normalized rates divided by the corresponding vari- 
ances is distributed as a x^ with a number of considered 
phase bins minus one degree of freedom. In all energy bands, 
we found that the normalized pulse profiles are equal within 
a statistical uncertainty of one standard deviation (see e.g., 
Fig.[l| where xL = 1-18 for 99 d.o.f.). We can estimate our 
sensitivity measuring a deviation from the template pulse 
profile, based on the assumption that the flux increase is 
constant throughout the phase. In the 20-40 keV energy 
range, we would be able to appreciate an increase in the 
un-pulsed emission as small as 4% at the 99% confidence 
level; our sensitivity diminishing to 6% in the 40-80 keV 
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Figure 1. Upper panel: pulse profiles measured by 
IBIS/ISGRI in the 20-40 keV energy band: crosses and di- 
amonds represent the pulse profile before and during the 
Crab flare, respectively. Lower panel: significance in stan- 
dard deviations of the difference between the pulse profile 
during and before the flare. 



band, and 10% between 80 and 150 keV. Above 150 keV, 
the sensitivity is limited to 30%. 

2.2. FERMI 

FERMI/LAT is a pair-conversion telescope that operates in 
the 30 MeV-300 GeV e nergy range with unprecedented sen- 



sitivity and resolution ( Atwood & et al. 2009 ) . In our analy- 
ses, we used all the events in the 0.1-300 GeV energy range 
belonging to the "diffuse" class, but rejected photons with 
zenith angles larger than 105° to avoid contamination from 
the Earth's bright 7-ray albedo. We limited our sample to 
the data taken between 2010 September 1 to October 6. The 
analysis was performed using the ScienceTools, provided by 
the Fermi collaboration (Version v9rl5p2). To obtain the 
source flux, we used the g tlike tool, which ex ploits a max- 
imum likehhood method ( [Mattox et"aI1|1996[) . The instru- 
ment response function is P 6 _V3 .DIFFUSE and the Galactic 
emission is reproduced using the model gll_iein_v02 • f itj^ 
All sources listed in the Fermi-LAT first year catalogue 
( [The Fermi-LAT Collaboration|2010] ) within 7° of the Crab 
position and with a TS detection greater than 10 are taken 
into account in the likelihood analysis. Each point in the 
source light-curve (Fig. [2| is obtained by means of a single 
likelihood analysis. Owing to the low number of events col- 
lected in every time bin, we used a PowerLaw^ to represent 
the total emission of the Crab. To optimize the signal-to- 
noise ratio (S/N), we accumulated the light-curves with dif- 
ferent binnings (6h, 12h, Id, 2d), finding that the 12-hour 
binning is the most adequate to follow the Crab count-rate 
evolution. In Fig. [2| the green dashed line indicates the av- 
erage fl-ux fr om the Crab of (286±2) x 10~^ ph cm~^ s~^, as 
reported by 
clearly visib 



Buehler et al. ( 2010[ ). Three distinct peaks are 
e, which exceed by a factor of nearly three the 



pre-flare flux level, corresponding to a detection significance 
> 15ct. The first steep increase occurs between 6:00 and 
18:00 UTC on September 18. In 12 hours, the flare reaches a 



http : //www. jb .man. ac .uk/~pulsar/ crab .html 



^ http://femii.gsfc.nasa.gov/ssc/data/analysis/software/ 
/aux/gllJem_v02.fit 

^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/ 
/xml_modeI_defs.html7^powerlaw2 
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Figure 2. Light curve of the Crab in the range 0.1-300 
GeV with 12h time bins. The green dashed hne indicates 
the flux of the Crab estimated over all the Fermi operation 
period. The green solid line is the average flux between 



Buehler et aPpOTO l. 
ux before and during 



September 19 and 21, as reported by 
The red lines represent the average f 

the flares (see Table [I]). The time is expressed in Modified 
Julian Date 



maximum flux of (8. 6±1.1) X 10 ^ ph cm ^s ^. The second 
flux increase occurs about two days after, between 12:00 
and 24:00 on September 20. The peak flux is (9.4 ± 1.0) x 
10~^ ph cm~^ s~^. The last flux increase was detected be- 
tween 18:00 on September 21 and 18:00 on September 22, 
reaching a value of (7.5 ± 1.1) x 10"^ ph cm~^ s~^. In all 
cases, the flare decay time was about 1 day. On the basis 
of the lightcurve reported in Fig. [2j we divide the data into 
three different sets as reported in Table [l] 



Table 1. Time intervals used to perform the Fermi analysis. 



Pre-flare 



Flare 



Post-flarc 



Start 
[UTC] 



2010-09-01 
00:00:00 



2010-09-18 
00:00:00 



2010-09-26 
00:00:00 



Stop 

[UTC] 



2010-09-16 
23:59:59 



2010-09-23 
23:59:59 



2010-10-05 
23:59:59 



In Fig. [Sj we compare the pulse profile measured before 
and during the flares. We can estimate the contribution of 
the nebula, by assuming that the pulsar emission is negli- 
gible in the phase ranges (j) < 0.25 and (j) > 0.8. Before the 
flares, the average rate of the nebular emission is (2.44 ± 
0.29) X 10~® ph cm~^ s~^ biUpj^g^^g, whereas during the flare 
it is (8.12 ± 0.55) x 10"* ph cm'^ s"^ bin-j^^^^^. After sub- 



tracting the contribution from the nebula, we compute the 
average rate of the pulsed emission: before the flare it is 
(9.32 ± 0.66) X 10"* ph cm'^ biUpj^^^^, whilst during 

the flares it is (8.52 ± 0.90) x 10"^ ph cm^^ g-i bin"!^^^^^. 
During the flares, we conclude that the flux from the nebula 
increased by a factor 3.33±0.46. In contrast, the pulsar flux 
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Figure 3. Exposure corrected Crab pulse profiles. Left 
panel: before the fiares. Right panel: during the flares. The 
horizontal green and red lines indicate the off-pulse aver- 
ages. 

remained constant within the errors. This result conflrms 



the preliminary claim of Hays et al. (20101. 

To study the high energy spectrum of the Crab, we ex- 
ploit the maximum lik elihood method for t he dataset listed 
in Table [I] Following |Abdo et all ( |2010[ ), we model the 
spectral emission from the Crab using two power-law com- 
ponents for the nebula, which represent respectively the 
synchrotron and IC emissions, and a power-law with an 
exponential cutoff to describe the pulsar contribution. As 
the flare is not pulsed, we fix the paramete r relat ed to 
the pulsar emission to these of Abdo et al. (2010), and 
assume that the IC contribution of the nebula has not 
changed during the fiare. The only free parameters are 
the synchrotron contribution and the normalization of the 
Galactic emission. The synchrotron emission is found to 
increase from a flux of (5.6 ± 1.3) x 10~^ph cm~^ to 



(32.4 ±2.7) X 10-^ ph cm 



in the 0.1-300 GeV band. 



3. Discussion 

To discriminate among the various models capable of repro- 
ducing the quasi-exponential turnover of the synchrotron 
emission of the nebula that peaks below the LAT en- 
ergy window, we studied the Fermi data complemented by 
archival CGRO/COMPTEL data (0.75-30 MeV). A single 
power-law cannot reproduce the pre-flare data (x^/d.o.f. ~ 
44/15). A power-law with a high energy exponential cut- 
off can instead reproduce the data {x^/d.o.f. ~ 4/14). To 
model the nebular synchrotron spectrum, we used the fol- 
lowing function 



dA^ 
dE 



— y 

IGcV/ 



exp 



E 



E, 



cutoff 



(1) 



The best-fit solution yields F = -2.20 ± 0.08 and 
A^o = (4.3±1.9)-10-i'' ph cm-2 s"i MeV-^. The difference 
between the quiescent and flaring spectra can be under- 
stood by considering two different extreme cases of either a 
constant power-law normalisation or a constant cutoff en- 
ergy. In the former case, an increase in the energy cutoff of 
a factor of nearly 5 (from 77 ± 15 MeV to 367 ±45 MeV) is 
needed (as illustrated in Fig. Ill) . This increase is averaged 
over the whole flaring period, thus represents a lower limit, 
since in each single flare the cutoff energy might have been 



3 



M. Balbo et al.: Twelve-hour spikes from the Crab Pevatron 



> 10 




E [GeV] 

Figure 4. Crab spectral energy distribution in the 100 MeV 
- 300 GeV energy range. The points with error bars are 
the Fermi detections before the flare (dark green), during 
the flare (red), and after the flare (light green). The black 
dashed line represents the contribution from the pulsar. The 
black dot-dot-dashed line represents the IC emission from 
the nebula. The blue and magenta dot-dashed and solid 
lines are the synchrotron nebula and the total emission be- 
fore and during the flare, respectively. Arrows indicate the 
95% confidence flux limits. 



higher. In the latter case, the spectral variability can be ex- 
plained by raising the continuum normalization by a factor 
of - 5. 

We note that the non-detection of any significant hard 
X-ray variability during the flare does not allow us to dif- 
ferentiate between the two possibilities as several electron 
po pulations are probably presen t in the nebula. As reported 
by [Aharonian fc Atoyan] ( |1998[ ), the COMPTEL data are 
characterized by a flattening of the spectrum that can be 
ascribed to the synchrotron emission of a separate electron 
population confined in compact regions such as wisps or 
knots. The luminosity of this component is less than 1% of 
that of the whole nebula. 

The cutoff in the synchrotron spectrum occurs at a char- 
acteristic frequency fpeak 4.2 x 10^ i? 7^. Provided that 
synchrotron radiation is the dominant mechanism through 
which particles channel their energy, the maximum elec- 



tron Lorentz factor obtained by equating tgync to taccci is 
7max oc (Brj)"-^/^, where 77 > 1 is the gyrofactor that char- 
acterizes the acceleration rate 7accoi = 7/^accci and iaccci — 
rjE/qeBc. This makes I'poak independent of B, leading to 
an electron synchrotr on energy cutoff « leOry-i MeV (see 
e.g. Aharonian] 2000 ) . A higher value may imply that the 
conditions in the accelerator differ from those in the emis- 
sion region, e.g. there is a lower magnetic field in the former, 
or that the synchrotron gamma-rays are produced in a rela- 
tivistically moving region, which produces a shift in the en- 
ergy cutoff to higher energies by the corresponding Doppler 
factor S. In this scenario, a value S — 367/160 « 2.3 would 
be required to explain the energy cutoff obtained during 
the flaring episode. 

On the other hand, magnetic fields at the level of be- 
tween ~ 300 fiG and — 2 mG are found in t he syn- 
wisps, respectively (see e.g 

1 GeV implies 



7 — 3— lOxlO^in the emitting regions. Taking 5 — 2.3, 
the comoving cooling timescale for those particles, taking 
an extreme value B ^ 2 mG, is — 0.3 d. The corresponding 
observer timescale would then be similar to the decay time 
of the peaks present in the Fermi lightcurve during the flar- 
ing period, < 1 d. In contrast, the flares could be related to 
an enhanced electron population, and the spectral variabil- 
ity could be obtained by raising the continuum normaliza- 
tion by a factor of ~ 5 or by adding a hard very high energy 
electron population (leading to a photon index F < 1). In 
this case, a Doppler boosting may not be required, and the 
observed duration of the flares could correspond to the syn- 
chrotron timescale of PeV electrons embedded in magnetic 
fields < 1 mG. 

The duration of the three short flares limits the size 
of the emitting region(s) to < lO'^^ cm. The peak lumi- 
nosity of these flares is higher/brighter than 10'^^ erg/s, 
i.e. ^ 0.5 %o of the Crab spin-down luminosity, assuming 
an isotropic distribution. The distance between the emit- 
ting region and the pulsar can thus be constrained to be 
^ 6 X 10^^ cm, i.e. not larger than 15% of the size of the 
bright synchrotron torus observed by Chandra and HST, 
and probably consistent with the half-width of this torus. 
The emitting region could therefore be linked to the inter- 
action zone between the jet and the torus, which is found to 
have brightened in th e HST image obtained on 2 October 



chrotron nebula and 



2008 1. Synchrotron radiation at 



Hester 
that 



(Caraveo et al. 2010). The three flares separated by two 
days could possibly be related to various emitting knots in 
this region. Alternatively, gamma-rays could be produced 
within the jet itself. However, if the emitting region were 
moving at relativistic speeds, the emission would be radi- 
ated within an angle — 1/(5. For reasonable values of the jet 
inclination angle wit h respect to the line of sight (see e.g. 
|Ng & Romani 20041, this scenario would make the flares 
difficult to detect. 

To conclude, the flare relative short durations (< 1 day), 
their soft spectrum, and the analysis of the pulse profile in 
the 0.1-300 GeV indicate that one or more compact por- 
tions (< 10^^ cm or < 0.1") of the synchrotron nebula are 
responsible for the flares. In these region(s), freshly acceler- 
ated PeV electrons are rapidly cooling, causing the observed 
variability. 
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